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Abstract 



Monte Carlo results for the low-lying glueball spectrum using an improved, 
anisotropic action are presented. Ten simulations at lattice spacings ranging 
from 0.2 to 0.4 fm and two different anisotropies have been performed in order 
demonstrate the advantages of using coarse, anisotropic lattices to calculate 
glueball masses. Our determinations of the tensor (2++) and pseudovector 
(1^ ) glueball masses are more accurate than previous Wilson action calcu- 
lations. 
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I. INTRODUCTION 



Numerical simulations of gluons on a space-time lattice provide at present the most 
reliable means of studying glueballs. Glueball correlation functions are, however, notoriously 
difficult quantities to measure in Monte Carlo simulations: since the masses of these states 
are rather high and their creation operators have large vacuum fluctuations, the signal-to- 
noise ratio falls extremely rapidly as the separation between the source and sink is increased. 
Because of this, reliable studies of glueballs on fine lattices have required prohibitively large 
computer resources. Thus, the development of more efficient simulation techniques in lattice 
QCD is crucial to establishing a detailed description of glueballs and their interactions. 

The objective of this work is to examine the effectiveness of using an improved, 
anisotropic lattice action to reduce the computational effort needed to determine the glue- 
ball spectrum in quenched QCD. Improved actions allow access to continuum physics on 
coarser lattices than possible using the simple Wilson discretization. Coarse lattice simula- 
tions are more efficient for several reasons: for a given physical volume, much fewer lattice 
sites are needed; the alleviation of critical slowing down permits the faster generation of sta- 
tistically independent gauge-field configurations; glueball operator smearing is faster due to 
the decreased number of links and a decrease in the number of smearing iterations required; 
glueball wavefunctions extend over much fewer lattice sites on a coarse lattice, making the 
variational technique far more effective when using a feasible number (a dozen or so) of basis 
operators. 

However, for glueball mass calculations, the coarseness of the temporal lattice spacing 
is a severe drawback. As the masses in lattice units of the states of interest are so large, 
the number of correlator time intervals which can be measured is reduced greatly A 
straightforward solution to this problem which preserves the computational advantages of 
coarse lattices [@-^ is to make use of anisotropic lattices in which the temporal spacing is 
much smaller than that in the spatial directions. This enables us to exploit the enhanced 
signal-to-noise of the correlation functions at smaller temporal separations. A natural scale 
for the temporal lattice spacing should be the inverse of the energy of the states of interest; 
thus, for glueballs, a temporal cut-off larger than 1.5 GeV allows resolution from accessible 
statistics of the correlator over a few time-slices. Meanwhile, the scale for the spatial lattice 
should be set by the size of the wavefunction of the state; a spatial grid separation in the 
range 0.2 — 0.4 fm would seem reasonable. 

Since we propose to use lattices in which the temporal lattice spacing is small, improve- 
ment of the discretization in this direction is not needed. Thus, a lattice action which couples 
only nearest-neighbor time-slices can be used. The transfer matrix corresponding to such an 
action is Hermitian and positive definite; all of our effective masses must converge to their 
plateau values monotonically from above. This ensures the validity of variational techniques 
which minimize the effective masses at small temporal separations. Such techniques are very 
effective in diminishing the excited-state contributions to the glueball correlation functions 
and are crucial for efficient extraction of ground-state masses. 

In this paper, we demonstrate the increased efficiency of glueball simulations using these 
actions on anisotropic lattices. We present results for the masses of three of the lighter 
SU(3) glueball states, the scalar (0+"*"), the tensor (2+"'"), and the pseudovector (l"*""). The 
masses of the first excited states in the scalar and tensor channels were also examined. 
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Ten simulations at lattice spacings ranging from 0.2 to 0.4 fm were performed, enabling 
reliable extrapolations to the continuum limit (although the mass of the scalar glueball was 
somewhat problematic). The results are compared to previous simulation data obtained 
using the Wilson action and we find that more accurate determinations of the tensor and 
pseudovector glueball masses have been achieved. A comparison of efficiencies is also made. 
Lastly, finite volume effects are shown to be small. 

The new action used in our simulations is described in Sec. y. The details of the glueball 
simulations, including the construction of the glueball operators, the generation of the gauge- 
field configurations, and the analysis of the Monte Carlo data, are given in Sec. The 
hadronic scale tq is used to relate our results at different values of the coupling j3 and 
the aspect ratio ^. The determination of this scale in terms of the lattice spacing using the 
static potential is outlined in Sec. |rV|. Sec. [V| contains our results and discussion: the glueball 
mass measurements are presented in detail; finite volume effects are studied; extrapolations 
of the masses at finite spacing to the continuum limit are undertaken; the conversion of our 
results into physical units is described; and a comparison of efficiencies with Wilson action 
simulations is made. Our conclusions are given in Sec. [V^, along with an outline of future 
work. 



II. AN IMPROVED, ANISOTROPIC DISCRETIZATION OF QCD 



Our glueball mass determinations rely on numerical simulations of glueballs on a Eu- 
clidean space-time lattice with spatial and temporal spacings and at, respectively. The 
improved gluonic action used in this study is given by 



S 



II 



P 



5 ^sp 4 ^^tp 1 ^sr 1 'C^str 

Sfwf 3v^t ~ 12 ^ ~ 12 uiu^t 



(1) 



where f3 = Q/g"^, g is the QCD coupling, Ug and Ut are mean link renormalization parameters, 
C, is the aspect ratio = ag/at at tree level in perturbation theory), and flc = J2c 1^-6 Tr(l — 
Wc), with Wc denoting the path-ordered product of link variables along a closed contour C 
on the lattice, flgp includes the sum over all spatial plaquettes on the lattice, fltp indicates 
the temporal plaquettes, Qsr denotes the product of link variables about planar 2x1 spatial 
rectangular loops, and figtr refers to the short temporal rectangles (one temporal link, two 
spatial). Explicitly, 



EElReTr 

X i>j 



- Ui{x)Uj{x+i)Uj{x+j)Uj{x) 
^tp = EE l^^eTr [l - Ut{x)U,{x+i)Ul{x + i)Ul{x) 

X i 

Qsr = EE l^^eTr [l - Ui{x)U,{x+i)U,{x + 2i)Uj{x + i+j)Uj{x+j)U^^{x) 
^str = EE |l^eTr [l - Ui{x)U,{x + i)Ut{x + 2i)Uj{x+i+i)Uj{x+i)U}{x) 



(2) 
(3) 

(4) 

(5) 



where x labels the sites of the lattice, i,j are spatial indices, and U^{x) is the parallel 
transport matrix in the gluon field from site x to x+jl. 
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This action, intended for use with at <^ ag, has 0{ I) discretization errors. The 

0{af) errors can be removed by the addition of counterterms which couple next-nearest- 
neighbor time-shces, but this introduces spurious high-energy modes which can cause consid- 
erable problems for our glueball simulations. These unphysical states appear in perturbation 
theory as additional poles in the gluon propagator. Their detrimental effects on the glueball 
correlation functions have been previously demonstrated Although these spurious states 
do not affect the asymptotic behavior of the glueball correlators, they do appreciably change 
the correlators at short temporal separations and can seriously hinder attempts to reduce 
excited-state contamination to hasten the onset of asymptotic behavior. Since our glueball 
mass measurements rely heavily on the reduction of such excited-state contributions to the 
glueball correlation functions, the use of an action which is free of spurious lattice modes 
is crucial. The action given in Eq. ^ couples only link variables on neighboring time-slices, 
which ensures that all of our effective masses converge to their plateau values monotonically 
from above and so validates the variational techniques employed. 

It is now known that perturbation theory by itself does not reliably determine the cou- 
plings in an improved action in lattice gauge theory. Hence, the interaction strengths in 
Sn have been determined using a judicious combination of perturbation theory and mean 
field theory. Mean field theory is introduced by separately renormalizing the spatial and 
temporal link variables: Uj{x) — > Uj{x)/us and Ut{x) Ut{x)/ut, where Ug and Ut denote 
the renormalization factors for the spatial and temporal links, respectively. The mean-link 
parameters Ut and Ug are best determined by guessing input values for use in the action, 
measuring the mean links in Landau gauge in a simulation, then readjusting the input values 
accordingly and tuning until the input values match the measured values. The determina- 
tion of these renormalization factors is described in more detail in Refs. 0,0. However, when 
at is significantly smaller than a^, we expect the mean temporal link Ut to be very close to 
unity since 1 — {^TrUt) oc {at/asY in perturbation theory. Hence, to simplify matters, we 
set Ut = 1. We introduce further simplifications by using a convenient and gauge-invariant 
definition for Ug in terms of the mean spatial plaquette given by Ug = {^KeTrPgg/)^^^, where 
Pggi denotes the spatial plaquette. This eliminates the need for gauge fixing, yields values 
for Ug which differ from those found using the Landau gauge definition by only a few per 
cent, and significantly speeds up the tuning process. 

At finite coupling g, the anisotropy ag/at is renormalized away from its input value ^. 
Measurements of this renormalization have been made using the static potential extracted 
from correlations along the different spatial and temporal axes of the lattice Without 
mean-link improvement, this renormalization can be as large as 30%. When the action 
includes mean-link corrections, this renormalization is found to be small, typically a few per 
cent. We used at/ag = ^ in all of our calculations, accepting the small radiative corrections 
to the anisotropy as finite lattice spacing errors, which vanish in the continuum limit. 

III. GLUEBALL SIMULATION DETAILS 

Glueballs may be labeled by their total (integral) spin J and their symmetries under 
spatial inversion and charge conjugation. However, on a cubic lattice, glueballs are charac- 
terized by their transformation properties under the cubic point group, combined with parity 
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and charge conjugation operations. The cubic group, Oh, has 24 elements that fall into five 
conjugacy classes, and thus, the dimensions of the irreducible representations (irreps) are 
1, 1, 2, 3 and 3. These irreps are labeled Ai, A2, E, Ti, and T2, respectively. Including 
parity and charge conjugation symmetry operations, there are 20 irreps (labeled by J^'" , 
where J now denotes an irrep of Oh)- In this study, four of the irreps which generate light 
(< 3 GeV) glueball states were simulated: the Af~^, E^^ , T-^^ , and channels. Of 
particular interest are the E~^^ and T^~^ irreps whose combined five rows correspond to the 
five polarization states of the tensor (2"^+) glueball which become degenerate as continuum 
rotational invariance is restored. This then gives information on the magnitude of lattice 
artifacts at finite cut-off. 

The mass of a glueball, G, having a given 

jPC 

can be extracted from the large-t behavior 
of a correlation function C{t) = (0|$^^''^(t) $'•■^^(0)10), where R denotes the lattice irrep 
corresponding to the J^'^ of interest and ^^^\t) = ^^^\t) — {0\^^^\t)\0) is a gauge- invariant, 
translationally- invariant, vacuum-subtracted operator capable of creating a glueball out of 
the QCD vacuum |0). As the temporal separation t becomes large, this correlator tends to 
a single decaying exponential limi^ooC*(^) = Zexp{—mGt), where mc is the mass of the 
lowest-lying glueball which can be created by the operator In order to extract mc, 

the correlator C{t) must be determined for t sufficiently large that C{t) is well approximated 
by its asymptotic form. However, the signal-to-noise ratio in any Monte Carlo determination 
of C{t) falls exponentially fast with respect to t. Thus, it is crucial to use a glueball 
operator for which C{t) attains its asymptotic form as quickly as possible. If \G) denotes 
the glueball state of interest, this means that we must choose an operator for which the 
overlap {G\^^^\t)\0)/[{G\G){0\¥^^^{t)¥^\t)\0)Y/^ is as near to unity as possible. For 
such an operator, the signal-to-noise ratio is also optimal [^. 

In order to construct such operators, we exploited the smearing and variational 
techniques which have been used with success in earlier Wilson action simulations. In 
each of the J^^ channels of interest, glueball operators were constructed on each time-slice 
in a sequence of three steps. First, smeared links Uj{x) and fuzzy superhnks Uj (x) were 
formed. Secondly, a set of basic operators (p^^ (t) were constructed using linear combinations 
of gauge-invariant, path-ordered products of the U^{x) and Uj (x) matrices about various 
closed spatial loops; each such linear combination was designed to be invariant under spa- 
tial translations and to transform irreducibly under the symmetry operations of the cubic 
point group according to the irrep of interest. Lastly, the glueball operators were 
formed from linear combinations of the basic operators, = Y^a'^^^'^'P^K't)^ where the 

coefficients v^J^^ were determined using the variational method. Each of these three steps is 
described below. 

Operators constructed out of smeared links and fuzzy superlinks have dramatically re- 
duced mixings with the high frequency modes of the theory. Thus, the use of spatially- 
smoothed links is an important part of reducing excited-state contamination in the glueball 
correlation functions. Two smoothing procedures were used: a single-link procedure and a 
double-link procedure. In the single-link procedure, every spatial link Uj{x) on the lattice 
is replaced by itself plus a sum of its four neighboring (spatial) staples, projected back into 
SU(3): 
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U;ix) = Vsu(3) {u,ix) + \s M^) U,{x + k) Ul{x+3)], (6) 

where Vsu{3) denotes the projection p|] into SU(3). Here, we denote this mapping of the 
spatial hnk matrices into the smeared hnk variables by sx^. In the double-link procedure, 
new superlinks of length 2as are built using neighboring staples which connect sites 
separated by a distance twice that of the length of the source link variables: 

U;{x)=Vsui3) [u,{x)U,{x + j) 

+ Xf Yl Uk{x) U,{x+k) U,{x+j+k) Ul{x+2j)\, (7) 

and we denote this mapping by fx^ Both procedures can be applied recursively; smeared 
links can be smeared again and fuzzy links of increasing length 2as, Aag, Sag . . . can be 
constructed. A smoothing scheme, S, is defined as a composition of single-link mappings 
and double-link mappings. Six different smoothing schemes were used. The simplest scheme 
used was the composition of two single-link smearings: Si = sx^ o sx^ ■ To simplify notation, 
we write this as Si = s\^ . We also used the compositions of four and six single-link mappings: 
1S2 = s\^ and 1S3 = s\^. In the other three smoothing schemes, the application of several 
single-link smearings, followed by one final iteration of double-link fuzzing was used: ^4 = 
fxf o ^5 = fxf o s\^, and Sq = fxj o s\^. Only one iteration of the fuzzing procedure 
which results in links connecting sites separated by 2as was found to be useful for the range 
of coarse Og values explored here. For the finer lattices (/3 = 2.4,^ = 5 and [3 = 2.6, C, = 
3), an extra four initial iterations of single-link smearing were used in all six smoothing 
schemes to enhance ground-state overlap. To simplify matters, the same values for the two 
parameters and A/ were used in all smearing and fuzzing iterations. These values were 
chosen to minimize excited-state contamination in the glueball correlation functions. A 
crude optimization was done in a set of low statistics runs and the optimal values A^ = 0.1 
and A/ = 0.5 were then used in all the glueball simulations. 

The second step in the construction of our glueball operators was the formation of a set of 
basic operators ^^^Ht) using linear combinations of gauge- invariant, path-ordered products 
of the Uj{x) and Uj (x) matrices about various closed spatial loops. Combinations which 
were Hermitian, invariant under spatial translations, and transformed irreducibly under the 
operations of the cubic point group according to the irrep of interest were constructed. For 
a more detailed exposition of this construction, see Ref. In each channel, a large set 
of prototypes were programmed, and a short simulation was then performed to determine 
the coefficients of each operator in the variational ground state. In each channel, the four 
operators with the highest of these contributions were then chosen for use in the production 
runs. The paths in this optimal set are illustrated in Fig. |l]. In the glueball simulations, 
these Wilson loops were measured on the link variables from the six smoothing schemes, 
yielding a total of = 24 basic operators 0^^^ (t) in each of the four channels. 

Finally, was formed from a linear combination of the basic operators, = 

Ea=ivi^^^^<^Kt)- The coefficients were determined using the variational method. First, 
the 24 X 24 correlation matrix was computed in the glueball simulations: 
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C'./.(t) = E(O|0i^nr+t)4''^(r)|O), (8) 

r 

where (l)^\t) denotes a vacuum- subtracted operator ^^f^t) = 0^^^^) - (0|</'i^n^)|0)- Note 
that (O|0^^^(t)|O) is independent of t. The coefficients w^^^ were then determined by mini- 
mizing the effective mass 



m{tj:,) = — In 



(9) 



where the time separation for optimization was fixed in all cases to to = 1- Let v'-'^^ denote a 
column vector whose elements are the optimal values of the coefficients . Then requiring 
dm{tD) / dv^^ = for all a yields an eigenvalue equation: 

C{tD) v(^) = e-™(*«) C{0) v(^). (10) 

The eigenvector Vq^^ corresponding to the largest eigenvalue e""^"*^*^^ then yields the coef- 
ficients VQa^ for the operator $Q^^(t) which best overlaps the lowest-lying glueball G in the 
channel of interest. Operators which overlap excited glueball states can also be constructed 
using the other eigenvectors of Eq. [10|. In particular, the operator ^\ (t) expected to best 



overlap the ffist-excited glueball state G* was obtained from the eigenvector corresponding 
to the second largest eigenvalue of Eq. 0. 

The elements of the correlator matrix given in Eq. ^ were estimated using the 
Monte Carlo method. Ten separate glueball simulations were performed on DEC Alpha- 
workstations. Configuration ensembles were generated using both Cabibbo-Marinari (CM) 
pseudo-heatbath and SU{2) sub-group over-relaxation (OR) methods. Link variables were 
updated in serial order on the lattice. We define a compound sweep as one CM updating 
sweep followed by three OR sweeps. In the glueball simulations, three compound sweeps 
were performed between measurements, and the measurements were averaged into bins of 
100 in order to reduce data storage requirements (except for the /3 = 2.6, ^ = 3 run in 
which 40 configurations were included in each bin). In all ten simulations, 100 bins were ob- 
tained. Our ensembles were tested for residual autocorrelations during the analysis phase by 
over-binning by factors of two and four; in all cases, the statistical error estimates remained 
unchanged. 

Values for the mean link parameter Us were determined self-consistently as previously 
described. This tuning procedure required a minimal amount of computational effort and 
provided thermalized configurations for later computations. The improved action simulation 
parameters used are given in Table |. 

For the data-fitting phase, the large 24 x 24 correlator matrices in each channel were 
reduced using the coefficients Vg^^ and v^'^^ to smaller 2x2 matrices CAsit) for A, B = 0,1: 

CABit) = E(0|5>f (r+t) ^f\rm. (11) 

The ground-state correlator Coo(i) was fit for t = tmin, • • • , ^max using a single exponential 

Coo(t) = Zoo {e-'"^* + e-'"^(^-*)} , (12) 
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where T was the temporal extent of the periodic lattice, to obtain an estimate of the mass 
me (in terms of a~[^) of the lowest-lying glueball in each channel. To determine the mass 
me* of the first-excited glueball and another estimate of mc, the 2x2 correlator was also 
fit for t = train, ■ ■ ■ , tma.x using the form 

CAB{t)= WBp{e-™-* + e-"^-(^-*)}. (13) 

p=G,G* 

Various fit regions tmin to tmax were used in order to check for consistency in the extracted 
values for the masses. Best fit values were obtained using the correlated method. Error 
estimates were calculated using a 1024— point bootstrap procedure; in all cases, error esti- 
mates were very close to being symmetric about the central best-fit values and were thus 
averaged to simplify presentation. 

IV. SETTING THE SCALE USING THE STATIC POTENTIAL 

In order to convert the glueball masses as measured in our simulations into physical units, 
we must set the scale by determining the lattice spacing at for each P and C, we consider. To do 
this, we must first choose one physical quantity to use as a reference. This reference quantity 
must then be measured on the lattice in terms of at- The experimentally- known value for 
the reference quantity is then used to extract the lattice spacing. A quantity which can 
be easily and accurately determined both experimentally and in numerical simulations is an 
ideal choice for such a reference. The mass of a low-lying particle is typically used for setting 
the scale. In our case, however, there are no unambiguous experimental determinations of 
the glueball masses, so instead, we must look for another purely gluonic quantity. 

The hadronic scale parameter tq defined in terms of the force between static quarks by 
[r^dV{r) /dr]r=ro = 1.65, where V{r) is the static-quark potential, is an attractive possibility. 
It can be measured very accurately on the lattice. The advantages in using tq to set the 
scale have been enumerated in Ref . |l^ . From phenomenological potential models, one finds 
tq ^ 0.5 fm. A disadvantage in using tq is that its physical value must be deduced indirectly 
from experiment, and there is some ambiguity in doing this, as will be discussed below. 
However, in the absence of a better gluonic reference, we have chosen rg to set the scale. In 
this section, we outline the determination of ro in terms of a^. 

In order to determine tq in terms of the lattice spacing, we need accurate measurements 
of the static-quark potential. We extracted V{f) for various spatial separations r, both on 
and off the axes of the lattice, from the expectation values of Wilson loops W{f,t) in the 
standard manner: 

W{r,t) = Z{r) ex\i[—tV{r)\ + excited state contributions. (14) 

In the Monte Carlo evaluation of the Wilson loops, measurements were taken after every 
four compound sweeps (as defined in Sec. The measurements of the Wilson loops were 
done independently of the glueball mass studies using separate ensembles of configurations. 
To minimize contamination from excited states, the Wilson loops were constructed from 
iteratively smeared spatial links. The single-link smearing method described previously 
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was used. A given smearing scheme is specified not only by the parameter A^, but also 
by the total number of smearing iterations, denoted by nx- Two different choices of the 
smearing parameter were used in all cases: one smearing was chosen to work well for small 
r = 1?^, the other to work well for large r. Separate measurements for each smearing were 
taken; cross correlations were not determined. The statistical noise in the evaluation of 
W{f,t) was reduced dramatically, especially for large temporal separations, by constructing 
the Wilson loops, whenever possible, from thermally-averaged temporal links [|r^. The 



thermal averaging was accomplished using the Cabibbo-Marinari pseudo-heatbath method 
(40 updates). Other relevant run parameters are given in Table 

The values of the potential V{r) were extracted from the Wilson loop measurements by 
fitting W{r, t) to the exponential form Z{f) exp[— tV^(r)] in the range t = tuim, • • • , ^max, for 
each f. The plateau region from tmin to tmax was chosen separately for each f in order to 
minimize the uncertainty in the extracted values for V^(r) while maintaining a good quality Q 
of fit. Best fit values were determined using the standard x-square test, taking into account 
temporal correlations among the W{f, t). The covariance matrix in was determined using 
the jack-knife procedure, and estimates for the uncertainties in the extracted values for V{r) 
were computed using the bootstrap method. Binning of the data was done as a crude check 
that our measurements were statistically independent. The results of a typical fit are shown 
in Fig. 1^, which is an effective mass plot for V{r) for f/ag = (2, 2, 2). The effective mass for 
V{r) is a function of t, defined as \ia[W{r,t)/W{r,t + at)], which tends to the true mass as 
t becomes large. 

Once a suitable plateau region in the effective mass was established for each V^(r*), the 
hadronic scale ro/a^ could be determined. We found that the on-axis potential V{r) for 
the range of r values studied here using coarse lattices fit a Coulomb plus linear form 
V{r) = Cc/r + err + Vq very well (with qualities of fit ranging from Q = 0.25 to Q = 0.99). 
We, therefore, used this form to interpolate V{r) and the force between static quarks. 
Simultaneous fits of the Wilson loops for the on-axis potential to the form Z{r) exp[— t(ec/r+ 
crr-f Vo)] were done, taking into account all correlations among the W{r,t) for both different 
t and r. Different regions in t were used for different r values; the plateau regions determined 
previously were used. Only the on-axis potentials were used; this prevented the covariance 
matrix in the to be minimized from getting too large. This covariance matrix was 
evaluated using the jack-knife method; uncertainties in the fit parameters Cc, cr, Vq, and 
Z{r) were obtained using the bootstrap method. Once we had an ensemble of bootstrap 
estimates for these fit parameters, the ratio ro/os and its bootstrap uncertainty were then 
determined using 



To/ a, = \J (1.65 + ec)/aa2. (15) 

Note that to compute ro/ag, we need the ratio as/ at since our fits yielded estimates of atV{f) 
only. We used the input value ^ since we know that its renormalization is small. Results for 
ro / as are given in Table |IT1[ 



Using the results in Table |T|, we can now express all energies measured in simulations 
in terms of tq. For example, in Fig. ^, we show the potential, including off-axis inter-quark 
separations, expressed in terms of tq. Lattice spacing errors are seen to be small. 
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V. RESULTS AND DISCUSSION 



A. Glueball mass measurements 

To allow clear resolution of the scaling properties of the low-lying glueball masses in the 
improved action, two sets of simulations were performed at two different anisotropies: six 
lattice spacings for an aspect ratio ^ = 3 and four spacings for = 5 were studied. The 
input parameters used in these simulations are given in Table |. 

The results of fitting the variationally-optimized correlators C{t) to the functions given in 
Eqs. |l^ and |13| are summarized in Tables |rV| - PCIIl| . Effective mass plots for the two smallest- 



ttg simulations are presented in Figs. and Figs. H— ^for ^ = 3 and ^ = 5, respectively. 
For each channel in each of the ten simulations, it was possible to find a fit region tmin ~ ^max 
in which the correlation function was well described by its asymptotic form as indicated by 
the quality of fit. In other words, convincing plateaux were observed in all effective masses. 
The most impressive plateau, observed in the Af^ channel for f3 = 2.4 and ^ = 5, spanned 
ten time-slices. In most cases, the onset of the plateau occurred when the source and sink 
operators were separated by only one time-step. The overlaps with the lowest-lying states 
were also found to be extremely good, better than 90% in most cases and often consistent 
with unity. This clearly demonstrates the effectiveness of the link-smearing and variational 
techniques in diminishing excited-state contamination. Fits using tmin = 2 or were also 
done to check for agreement with the tmin = 1 results. At time separations for which the 
ground state could be reliably observed, the off-diagonal elements of the reduced correlation 
matrices CAsit) were found to be consistent with zero within statistical uncertainty. This 
suggests that the link-smearing, variational method also gives an excellent construction of 
the first-excited state in each channel. 

Our best estimates for the glueball masses in terms of a^^ are indicated in boldface in 
each of the Tables [rV[ - pCllI| . These estimates are summarized in Table pClV| . Masses for the 



first-excited states are also indicated in the Ng^p = 2 fits listed in the Tables fVl - pCIII 



B. Finite volume effects 

In this work, we were concerned with the magnitude of discretization errors in the glueball 
mass determinations from coarse lattice simulations using an improved action. In order to 
evaluate these errors, we had to eliminate uncertainties from all other sources. The increased 
efficiency of simulations on coarse, anisotropic lattices allowed us to reduce statistical errors 
to the acceptable level of about 1%. The only remaining source of uncertainty we had to 
address was finite volume. The masses of particles confined in a small box with periodic 
boundary conditions can differ appreciably from their infinite volume values; finite volume 
effects can also induce a splitting in the masses of the E and T2 tensor polarizations. Finite 
volume effects on the scalar glueball mass have been analyzed before []12| , but the effects on 
the tensor and the pseudovector are less well known. 

In order to ascertain the effects on our glueball masses of simulating in a finite volume, 
four extra simulations were performed for j3 = 2.4, ^ = 3 using lattices of spatial extent 
Ls/tts = 6,5,4, and 3. The temporal extent was held fixed at 24 grid points. For each of 
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these volumes, the mean-field renormahzation parameter Us was recalculated. The 3^ lattice 
was the only simulation that required any change in this parameter, and in this case, the 
effect was small; Ug increased by only 0.3%. The results from the L^/as = 8, 6, 5, and 4 runs 
for the glueball masses in terms of aj~^ are given in Table PCV| . Note that the results from the 



6'^ lattice differ very little from those from the 8^ lattice, suggesting that our lattice volumes 
are sufficiently large to ensure that finite volume errors are negligible. For the 3^ lattice, 
no plateaux in the effective masses for the Af~^, E^^ and Tg*"*" channels were observed; the 
mass in the Ti~ channel was found to be 1.44(1). The operators used in these runs were 
the same as those constructed for the large volume runs and thus, were not optimized to 
give large overlaps with the light torelon states present in small volumes. It is likely that 
this effect was responsible for the poor overlap of our operators with the scalar and tensor 
states on the 3^ lattice. 

The properties, such as the mass, of a glueball confined in a small box with periodic 
boundary conditions differ from those in an infinite volume. The modification of the mass 



of a particle due to finite volume effects has been estimated in Ref. |]13 
atmciz) = atmcioo 



1 — Ag exp {-V^z/2)/z\ , (16) 



where z is the dimensionless length scale z = 171^++ Lg, m^++ is the infinite- volume mass 
of the scalar glueball, and Ag is related to the strength of an effective triple scalar glueball 
interaction vertex. The mass shift given in Eq. |16| is valid for sufficiently large z and arises 
from the exchange of scalar glueballs across the periodic boundaries of the lattice. Finite 
volume errors in our glueball masses measured on an 8^ lattice at /3 = 2.4, ^ = 3 (where Lg fti 
2 fm, similar to the volumes used in the other nine simulations) can be estimated by fitting 
the form given in Eq. |16| to the masses in Table PCV| . Let u = atm^++ and p = ^Ls/ag, then 

z = poj. The fit was done first using the function atm(z) = to — \^++ exp(— -\/3pco'/2)/p, 
where u and A^++ were the fitting parameters. The best fit value for u was then used in 
the fits to the results for the other irreps; to simplify matters, the uncertainty in u was 
neglected in these fits. The fit also included the energy estimate extracted from the 
Ls/tts = 3 simulation. The results of these fits are summarized in Table pCVl|; the estimates 



of the finite volume errors are listed in the final column of this table and are given by 
maiSuj^) /mG{oo) — 1 using Eq. 0. In all cases, these errors were insignificant compared 
to statistical errors; this means, for example, that any differences between the large-volume 
masses in the T2 and E channels must be due purely to discretization errors. It is interesting 



to note that our estimate of A^++ agrees well with the value 190 ± 70 found in Ref. |T2 



C. Continuum limit extrapolations 

The glueball mass estimates in terms of a^^ were combined with the determinations of 
the hadronic scale ro/ag. The results are shown in Figs. |T2|, |13|, and |l^. In these figures, 
the dimensionless product of ro and the glueball mass estimates are shown as functions of 
(a^/ro)^. Solid symbols indicate results from the ^ = 3 simulations, while open symbols 
are used for the results from the ^ = 5 runs. In Fig. 12, the lowest-lying masses in each 
of the channels Ai~^, E'^^, ^2^^; and Ti are compared with results from small-Og Wilson 
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action simulations. The lowest-lying and first-excited masses in the E~^^ and T^^ channels 
are shown in Fig. 13, and the ground-state and first-excited state in the Af^ channel are 



depicted in Fig. [T^. To extract physical predictions (for the pure-gauge theory), the curves 
in these plots must be extrapolated to the continuum limit a^/rQ 0. Discretization errors 
are given by the deviations of the finite-Os results from these limiting values. 

The lowest-lying states in the E~^^ and T-^^ channels correspond to the five polarizations 
of the tensor 2~^~^ glueball in the continuum. Differences between the E~^^ and T-^^ masses 
are a measure of violations of rotational symmetry due to finite spacing artifacts. In Fig. |1^, 
such violations are seen to be small for our less coarse lattices and become appreciable as 
the spacing gets very large. Discretization errors in the T2~^ exceed those of the E~^~^] on 
our coarsest lattices, finite spacing errors are only a few per cent for the E'^^ channel, but 
about 15% in the T^~^ channel. In the E~^^ channel, the ^ = 3 results differ very little from 
those using the higher aspect ratio ^ = 5, suggesting that the 0{af) errors are negligible. 
However, small differences between the results from the two anisotropics are visible in the 
Tj*"^ channel. One expects that 0{af) errors will decrease as ^ is increased. Since the Tg*"*" 
discretization errors are slightly larger for the ^ = 5 runs, 0{af) errors can account for this 
difference only if such errors offset the O(a^) errors. 

The leading discretization errors in the tensor glueball masses are expected to be 
0(a^, a^, a^a^). However, we have already argued that the results in Fig. ^ imply that 
the 0{af) errors are negligible. Since the action included mean-field correction factors, we 
also expected that 0(af) errors would dominate over 0(a<jO^) errors and in our continuum 
limit extrapolations, we assumed that this was true unless the fit provided compelling evi- 
dence to the contrary. Although we expected the leading discretization errors to be 0{a'^), 
the following three functions were used in our continuum limit extrapolations: 

V5o(a.) = roTJiG, (17) 
^2(0.3) = romc + C2 (a^/ro)^ (18) 
ipi{as) = roiriG + C4 (a^Ao)^, (19) 

where C2, C4 and romc are best-fit parameters. The results of these fits are given in Ta- 
ble PCVII| for the ^ = 3 data and Table PCVIII] for the ^ = 5 simulations. Comparing the 



values of x^/dof, one sees that the fitting function ip^ was preferred for both the ^ = 3 and 
^ = 5 results, although only marginally so for ^ = 3. Given this fact and our expectation 
that should best describe the leading discretization effects, we took rQinc from the (y94 fits 
as our continuum limit estimates (indicated in boldface in Tables pCVIl] and PCVIII|) . These 



four estimates are in very good agreement not only with one another, but also with the Wil- 



son action estimates. These fits using (p^ are shown in Fig. |T2[ For our final estimate of the 
tensor glueball mass, we performed a simultaneous fit with the four data sets (two irreps and 
two anisotropics) using four separate ^94 functions but constraining the intercept parameter 
rQ-rriG to be the same for all four fitting functions. This yielded rom(2"'"+) = 5.85 ± 0.02 
with x^/dof = 1.01, in agreement with the Wilson action estimate rQm{2~^~^) = 6.0 ± 0.1, 
obtained by fitting all of the Wilson action measurements shown in Fig. ^ to (pp. 

We also examined the discretization errors in the masses of the first-excited glueball 



states in the E~^~^ and T2'"'' channels. These are shown in Fig. [13. There are several 



reasons for interpreting these data as different polarizations of a spin-two excited state: 
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the two irreps extrapolate to the same continuum hmit value; if the Ta*"*" state were spin 
three, then there would be a degeneracy with the T^~^ and ^2"^ channels and this was not 
observed (these results will be presented elsewhere); if they were polarizations of a spin- four 
state, then again, a similar level must also be found in the Ti"^ channel. The degeneracy 
between the two irreps and the weak finite-volume dependence of their energies also makes 
an interpretation of this state as a torelon pair or a two-scalar-glueball scattering state 
unlikely, although the mass of this level lies close to twice the mass of the scalar glueball. 
Continuum limit extrapolations were performed using the three functions of Eqs. p!^^ - |T9| ; 
measurements from the two largest spacings for ^ = 3 and the single largest spacing for 
= 5 were excluded from these fits. The results are given in Tables pCVll| and pCVlll[ . Again, 
we expected 994 to provide the most reliable extrapolations to the — > hmit; this was 
confirmed by the fact that (f^ yielded E~^~^ and T^^ continuum limits in best agreement 
with each other. Differences found between the ^ = 3 and ^ = 5 extrapolations were not 
statistically significant. The fits using if 4 are shown in Fig. |T^. Our final determination of 
the mass of the excited tensor glueball, obtained from a constrained set of four (p^, similar to 
the ground-state extraction, was rom(2*++) = 8.11 ± 0.04, where x^/dof = 2.3. This mass 
has not been reliably determined in any previous simulations. 

Finite spacing errors in the mass of the pseudovector glueball were also studied. 
These were found to be small and are shown in Fig. |12|. The results from the different 
anisotropics are in good agreement. Extrapolations to the — ^ limit were done using the 
three functions ipo, ip2, and v?4; the results of these fits are summarized in Tables pCVIl] and 



XVIII| . The continuum limits obtained from fits to the ^ = 3 and ^ = 5 data agreed only for 
the constant fit form ipo. The fits to (p2 and yielded slope parameters (c2 and C4) with 
large relative errors and opposite signs for the different anisotropics. Hence, the function ipQ 
was used to extrapolate to the continuum limit. Due to the very good agreement between 
the ^ = 3 and ^ = 5 results, all ten data points were used in our extrapolation fit. Our 
estimate from this fit (shown in Fig. ^) was rom{l~^~) = 7.21 ± 0.02, where x^/dof = 1.55, 
in agreement with the extrapolation rom{l~^~) = 7.5 ± 0.4 using ipo of the Wilson action 
results shown in Fig. |l^. 

In contrast to the tensor and pseudovector, the scalar glueball mass showed significant 
finite spacing errors (see Fig. |12D, even for our less coarse lattices. As was increased, the 
scalar mass first decreased, reached a minimum near a^/ro ~ 0.6, then gradually increased. 
Near the minimum, the mass was about 25% lower than estimates of the continuum limit 
from small-Os Wilson action computations; a 20% discretization error was observed in the 
result from our smallest simulation. Although the magnitudes of these errors were sig- 
nificant, they were smaller than those obtained using the Wilson action by a factor of two. 
In order to extrapolate to the continuum limit, an appropriate fitting function was needed. 
The leading discretization errors were expected to be 0{af,a'^,asa'^). However, there were 
no distinguishable differences between the ^ = 3 and ^ = 5 results, suggesting that the 0{af) 
errors were negligible, leaving us to consider 0(a^, a^a^) effects. By inspection, one sees that 
the fitting form ip^, which neglects one-loop ©(a^a^) effects, cannot describe the data, in 
contrast to the data for the tensor and pseudovector glueballs. As — 0, we expect the 
coupling as{as) to vanish as —1/ In(asA), where A is an appropriate scale parameter. Hence, 
we were led to consider the following four-parameter fitting function: 
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(a IrnY 

VidcLs) = roTUG + C2 T^TT — -— ^ + C4 (a^/ro)^. (20) 

Cl - In [{as/roy\ 

However, it was not known how reliably the leading perturbative behavior of as{as) would 
describe the true cut-off dependence of the coupling over the large range of spacings consid- 
ered. Taking this into account and inspecting the behavior of the actual data, we decided 
to also consider the following simpler quadratic form: 

<y52,4(a.) = roniG + C2 {as/r^f + C4 {aJr^Y. (21) 

Both of these functions were fit to the mass measurements from all ten simulations; the 



results of these fits are summarized in Table [XIX| . The function (piL yielded a slightly better 



fit and a continuum limit for the scalar glueball mass of 3.98 ± 0.15. This fit is shown 
in Fig. |12|. An extrapolation of existing Wilson action data using yielded 4.33 ± 0.05. 
Given the quality of the scalar glueball mass estimates using the Wilson action, this slight 
discrepancy raises doubts concerning the reliability of the extrapolation using ifiL'-, mass 
estimates using the improved action for a few values of smaller than those considered 
here would be needed to resolve this discrepancy. 

One explanation for the 20% discretization errors in the scalar glueball mass is that the 
scalar glueball is extremely small. However, there is evidence [0,0 that the presence of a 



critical endpoint of a line of phase transitions (not corresponding to any physical transition 
found in QCD) in the fundamental- adjoint coupling plane is responsible for lowering the 
scalar glueball mass near the crossover region in the Wilson action. It is possible that the 
scalar glueball mass in the improved action used here may be similarly influenced. If so, 
the fact that this effect appears to be less pronounced for this action suggests the possible 
existence of other perturbatively-improved actions in which the scalar glueball mass is even 
less affected by scaling violations. We are currently searching for such actions. 

Discretization errors in the mass of the first-excited state in the A^^ channel were also 



found to be significant, as shown in Fig. 0. The mass of this state is nearly twice that of the 
lowest-lying scalar glueball, suggesting that this state may simply be two glueballs. Given 
the significant discretization errors in the single glueball scaling data, one would expect 
similar systematic errors in the two glueball state. The absence of any level of similar 
mass in all other channels justifies the spin-zero interpretation of this state. Considering 
the difficulties encountered in extrapolating the lowest-lying scalar to the continuum limit, 
we made no serious attempt to determine the continuum limit of this first-excited state. 
However, the result of a fit using Lp2,i is included in Table pCIX| . The possible interpretation 
of this level as a two glueball system might be strengthened by a more precise finite- volume 
study. 



D. Conversion to physical units 

In order to convert our glueball mass computations into physical units, we must specify 
the value of the hadronic scale. The hadronic scale tq has a precise definition in terms of 
the static-quark potential. However, the static-quark potential cannot be directly measured 
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in an experiment; it must be deduced indirectly from other observables. We decided to use 
a variety of different physical quantities to deduce tq. 

In Table \)iX[ estimates of r^^ using the results from various quenched lattice simulations 
are shown. For each computation, the quantity used to set the lattice spacing, such as the 
mass of the p or the IP — 15 splitting in heavy quarkonia, is indicated. The determination of 
Tq'^ from was accomplished using values of a/ro given in Ref. |jT6| for the Wilson gluonic 
action at various values of (3, interpolating where necessary. Note that due to quenching 
effects, To varies with the quantity used to set the scale. The entries in the last column of 
Table pCX] cannot be considered as different measurements of a single quantity and thus, 
strictly speaking, their weighted average has no statistical meaning. The last column of the 
table is meant to illustrate the range of values one obtains for r^^ when using various scale 
setting quantities. We expect that the value of r^^ appropriate for the low-lying glueballs 
should lie somewhere within this range. The simple average = 410 ± 20 MeV of the 
determinations in Table XX was taken as our estimate of the hadronic scale. 

For our final continuum mass estimates of the tensor glueball, we found 2400 ± 10 ± 
120 MeV (where the first error is statistical and the second is from uncertainties in the 
determination of tq). It is interesting to note that our mass estimate lies within 8% of the 
mass of the /j(2220) resonance |T^,|TB|, reported to have quantum numbers (even)"'""'". In 
order to make a direct comparison with experiment, however, corrections to our result from 
light quark effects and mixings with nearby conventional mesons must be taken into account. 
Our estimate of the mass of the first-excited glueball in the tensor channel was 3320±20±160 
MeV; for the pseudovector state, we found a mass of 2960 ± 10 ± 140 MeV. Our estimate 
from the fit using ipn for the mass of the scalar glueball was 1630 ± 60 ± 80 MeV; however, 
we regard the continuum limit extrapolation for this state as being less reliable than those 
for the other glueballs. 



E. Comparison of efficiencies 

A quantitative comparison of efficiencies is difficult to make. There are many factors 
which affect the overall efficiency of a Monte Carlo simulation. Certainly, the number of 
link updates is an important factor. The speed of an update is, of course, platform-, action-, 
and algorithm-dependent. On the DEC Alpha-workstations we used, a CM update using an 
improved action required twice as much time as for the Wilson action; the improved-action 
OR updating time was three times longer. Critical slowing down and thermalization are also 
contributing factors, but a crucial issue is the reduction of excited-state contaminations in 
the glueball correlators. In our coarse lattice simulations, we found that current methods for 
constructing good glueball operators were very effective in hastening the onset of plateaux 
in the effective masses. 

Given these difficulties in assessing the efficiency of a glueball simulation, we decided to 
make our comparisons based simply on the number of link updates Ni^ and the fractional 
error e attained in the final mass estimates. Since the error in a Monte Carlo estimate 
decreases with the number of measurements as 1/ ^/N, we expect that the reciprocal 
product of the number of link updates and the square of the fractional error is approximately 
proportional to the efficiency of a simulation; we denote this quantity hj S = l/{e^Niu). An 
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interesting comparison to make is between simulations at a small lattice spacing, such as 
Qs ~ 0.05 fm, using the Wilson action and improved-action simulations at a spacing ~ 0.2 
fm. Such a comparison is relevant because the discretization errors in these cases are of 
comparable magnitude, excepting the scalar glueball mass. For the (3 = 6.4 Wilson action 
run in Ref. |T9[ by the GFll collaboration using gauge invariant glueball operators, a total of 
3.13 X 10^2 link updates were performed, an error of 2.5% was achieved in the scalar glueball 
mass, and a fractional error of 3.6% was obtained for the tensor mass. Using tq = 0.48 fm 
and a^/ro = 0.101(2), the lattice volume in this simulation was (1.55 fm)^ x (1.45 fm). For 
the same value of (3 in Ref. |]2D| by the UKQCD collaboration, 1.35 x 10^^ link updates were 
made and fractional errors 3.4%, 3.3%, and 9% were obtained for the scalar, tensor, and 
pseudovector masses, respectively. The lattice volume was (1.55 fm)^ for this simulation. In 
our j3 = 2.6, ^ = 3 run, 5.76 x 10^ link updates were performed and 1.5%, 1.0%, and 1.0% 
errors in the scalar, tensor, and pseudovector masses were achieved. Our lattice volume was 
(1.93 fm)^ since a^/ro = 0.4021(9). For our f3 = 2.4, ^ = 5 run, the number of link updates 
was 9.83 X 10^ and the errors obtained in the scalar, tensor, and pseudovector masses were 
1.0%, 0.5%, and 0.8%, respectively. Using a^/ro = 0.459(1), our lattice volume for this run 
was (1.76 fm)^. Thus, the ratios of the S values for our /3 = 2.6 and P = 2.4 simulations to 
those of the GFll run were 1500 and 2000, respectively, for the scalar glueball mass, and 
7000 and 17000 for the tensor glueball mass. The ratios of the j3 = 2.6 and (3 = 2.4 £^-values 
to those of the UKQCD run were 120 and 160, respectively, for the scalar mass, 260 and 
600 for the tensor, and 1900 and 1700 for the pseudovector. 

Considering our ten simulations together, a total of 5 x lO^'' link updates were per- 
formed. For the Wilson action simulations of Refs. [|T9lE0i, an estimated 10^^ and 10^^ link 



updates were required, respectively, to generate continuum limit results whose statistical 
uncertainties were larger (for the tensor and pseudovector states) than those quoted here: 
the statistical error on our estimate for rom(2"*""'') was about five times smaller than that 
from the extrapolation of the Wilson action results, and for rom{l~^~), the uncertainty was 
twenty times smaller, implying that about 25 — 400 times greater statistics would be required 
for similar accuracy. Thus, in total, the anisotropic lattice simulations were certainly more 
than 1000 times more efficient. 

The above discussion illustrates the computational advantages of extracting non-scalar 
glueball masses from simulations on coarse, anisotropic lattices (a^ ~ 0.2 fm) using an 
improved action instead of lattices for which ~ 0.05 fm. The excellent overlaps achieved 
from our variational calculations demonstrate another advantage of simulating on coarse 
spatial lattices. The glueball wavefunctions extend over only a few points of the lattice 
when as ~ 0.2 fm. Thus, variational calculations using a feasible number of basis functions 
(a dozen or so) can yield very good approximations to the glueball wavefunctions. This 
will be especially important for future decay constant calculations and determinations of 
mixings with non-glueball states. These advantages have already enabled us to study the 
more massive glueball states which have yet to be simulated reliably using the Wilson action; 
these results will be reported elsewhere. 

These advantages are less clear for the scalar glueball mass due to the presence of 20% 
discretization errors at ~ 0.2 fm. Using the action of Eq. 0, simulations at one or more 
lattice spacings smaller than 0.2 fm would be needed to firmly establish the continuum limit. 
A more attractive approach would be to use an action for which discretization errors in the 
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scalar glueball mass at ~ 0.2 fm are negligibly small. The search for such an action is 
currently underway. 

VI. CONCLUSION 

We have demonstrated the advantages of using anisotropic lattices and an improved glu- 
onic action for simulating glueballs. Ten simulations at lattice spacings ranging from 0.2 to 
0.4 fm were performed using DEC Alpha-workstations, and the results were extrapolated 
to the continuum limit. Results for the masses of the scalar (O^"*"), the tensor (2++), and 
the pseudovcctor (1^ ) glueballs in SU{3) pure-gauge theory were presented in terms of the 
hadronic scale Tq^. The continuum limits for the tensor and pseudovcctor glueball masses 
were obtained with uncertainties of less than 1%, significantly improving upon previous esti- 
mates from Wilson action simulations carried out with the aid of super-computer resources. 
Extrapolation of the scalar glueball mass to the continuum limit was hampered by uncer- 
tainties in choosing the fitting function and discretization errors which were 20 — 25% even 
for our smallest lattice spacings; although uncomfortably large, these finite-spacing errors 
were half as large as those obtained using the Wilson action. Finite volume errors in our 
results were shown to be negligible. The masses of the first excited states in the scalar and 
tensor channels were also examined. 

Our results show that spatially-coarse, anisotropic lattice simulations are an effective 
means of studying gluonic systems. The techniques exploited here are sufficiently powerful to 
overcome the difficulties which plague Monte Carlo calculations involving gluonic excitations. 
These methods should be useful for studying the spectrum of heavier glueball states. Data 
for the masses of all twenty lattice irreps of the cubic group (including parity and charge 
conjugation) are currently being accumulated in order to survey the spectrum of SU{3) 
glueball states below 4 GeV comprehensively. We shall report on these results in the near 
future. We also plan to use the techniques outlined in this paper to determine various 
glueball matrix elements and decay strengths, to investigate the mixings of glueballs with 
conventional hadronic states, and to study mesonic states containing excited glue (the so- 
called hybrid mesons). The size of discretization errors in the scalar glueball mass was the 
only disappointing aspect of this work; we are currently investigating a new class of lattice 
actions with the hope of reducing these lattice artifacts for as ~ 0.2 fm to the level of a few 
per cent. 
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FIGURES 




FIG. 1. The four Wilson loop shapes in each channel used to form the lattice glueball operators. 
The complete set of 24 operators was formed by computing linear combinations of each of these 
loops rotated and translated across the lattice on six different sets of smoothed links. Where a 
loop shape occurs twice, it is used in two different projections into the appropriate irreducible 
representation. 
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FIG. 2. Effective mass plot sfiowing the results of a single-exponential fit to the Wilson loop 
for V{r) with r/ag = (2, 2,2), (3 = 2.4, and C = 3. The tmin — imax region of the fit is also indicated. 
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FIG. 3. The static-quark potential V{r) expressed in terms of the hadronic scale ro- This plot 
includes measurements from the /3 = 2.2, 2.4, and 2.6 simulations for ^ = 3, and the (3 = 2.2 and 
2.4 simulations for ^ = 5. Lattice spacing errors are seen to be small. 
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FIG. 4. Effective mass plot showing the results of a single-exponential fit to the glueball 
correlation function for the Af~^ channel for (3 = 2.6 and ^ = S. The tmin — ^max region of the fit 
is also indicated. 
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FIG. 5. Effective mass plot showing the results of a single-exponential fit to the glueball 
correlation function for the channel for P = 2.6 and ^ = 3. 
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FIG. 6. Effective mass plot showing the results of a single-exponential fit to the glueball 
correlation function for the E~^~^ channel for /3 = 2.6 and C = 3. 
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FIG. 7. Effective mass plot showing the results of a single-exponential fit to the glueball 
correlation function for the T2"'' channel for P = 2.6 and ^ = 3. 
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FIG. 8. Effective mass plot showing the results of a two-exponential fit to the 2x2 matrix of 
glueball correlation functions for the ^2^^ channel for /3 = 2.4 and ^ = 5. The imin — ^max region 
of the fit is also indicated. 
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FIG. 9. Effective mass plot showing the results of a single-exponential fit to the glueball 
correlation function for the channel for /3 = 2.4 and ^ = 5. 
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FIG. 10. Effective mass plot showing the results of a two-exponential fit to the 2x2 matrix 
of glueball correlation functions for the E~^^ channel for /? = 2.4 and ^ = 5. The tmin — ^max region 
of the fit is also indicated. 
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FIG. 11. Effective mass plot showing the results of a two-exponential fit to the 2x2 matrix 
of glueball correlation functions for the channel for /3 = 2.4 and ^ = 5. The tmin — imax region 
of the fit is also indicated. 
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FIG. 12. Glueball mass estimates in terms of rg against the lattice spacing (os/ro)^. Results 
from the ^ = 5 simulations for the lattice irreps Af^, , ^2"^ -^1 labeled o, □, O, and 

A, respectively. The corresponding solid symbols indicate the results from the ^ = 3 simulations. 
Data from Wilson action simulations taken from Refs. 



|]19|,^j2C,23| are shown using crosses. The 
dashed, dotted, and dash-dotted curves indicate extrapolations to the continuum limit obtained 
by fitting to the = 3 data, the ^ = 5 data, and all data, respectively. The solid line indicates the 
extrapolation of the Wilson action data to the continuum limit. 



30 



^ I I I I I I I I 

0.2 0.4 0.6 0.8 



FIG. 13. Estimates of the masses of the tensor gluebaU and its first-excited state in terms of 
ro against the lattice spacing (os/ro)^. Results from the ^ = 5 simulations for the E~^^ and T2~^ 
irreps are labeled by □ and O, respectively. The corresponding solid symbols show the results from 
the ^ = 3 simulations. The dashed and dotted curves indicate extrapolations to the continuum 
limit obtained by fitting to the = 3 and the ^ = 5, respectively (see Tables XVII and XVIII). 
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FIG. 14. Masses of the scalar glucball and its first-cxcitcd state in terms of tq against the lattiee 
spacing (os/ro)^. Open and solid symbols indicate results from the ^ = 5 and ^ = 3 simulations, 
respectively. The extrapolation to the continuum limit is indicated by the dash-dotted curve. 
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TABLES 



TABLE I. InpTit parameters Tised in the ghieball simulations. 



p 






Lattice 




4 


1.7 




3 


6^ X 18 




0.3075 


1.9 




3 


6^ X 18 




0.340 


2.0 




3 


8^ X 24 




0.356 


2.2 




3 


8^ X 24 




0.3885 


2.4 




3 


8^ X 24 




0.421 


2.6 




3 


10^ X 30 




0.45U5 


1.7 




5 


6^ X 30 




0.295 


1.9 




5 


6^ X 30 




0.328 


2.2 




5 


8^ X 40 




0.378 


2.4 




5 


8^ X 40 






TABLE II. Various 


run parameters for the static-quarlc potential measurements. 


including 


lattice sizes, 


total numbers of configurations used, and parameters for the two different 


smearing 


schemes. 












P 




Lattice 


# configs 


(Ai,nAi) 


(A2,nA2) 


1.7 


3 


6^ X 18 


2275 


(0.05,4) 


(0.10,4) 


1.9 


3 


6^ X 18 


1280 


(0.07,4) 


(0.14,4) 


2.0 


3 


8^ X 24 


848 


(0.07,4) 


(0.15,4) 


2.2 


3 


8^ X 24 


1024 


(0.10,6) 


(0.16,6) 


2.4 


3 


8^ X 24 


1024 


(0.12,6) 


(0.25, 6) 


2.6 


3 


10^ X 30 


1100 


(0.25,6) 


(0.50,6) 


1.7 


5 


8^ X 40 


810 


(0.10,6) 


(0.25,6) 


1.9 


5 


6^ X 30 


1024 


(0.08,4) 


(0.16,4) 


2.2 


5 


12^ X 48 


315 


(0.20,4) 


(0.30,4) 


2.4 


5 


8^ X 40 


548 


(0.20,6) 


(0.40, 6) 
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TABLE III. Results for the hadronic scale tq in terms of the lattice spacing Og. The Coulombic 
coupling e,. and the string tension a obtained from a fit of the on-axis potential to a Coulomb plus 
linear form V{r) = Cc/r + ar + Vq are also given. 



(3 


i 


ro/as 


as/ro 






2 


1.7 


3 


1.161(2) 


0.861(2) 




-0.07(1) 


1.58(1) 


1 Q 










—0 1 '^(9) 


1 59('9'l 


2.0 


3 


1.375(1) 


0.7271(8) 




-0.188(7) 


1.462(7) 


2.2 


3 


1.615(2) 


0.6192(8) 




-0.288(8) 


1.362(8) 


2.4 


3 


1.978(6) 


0.505(1) 




-0.321(6) 


1.329(6) 


2.6 


3 


2.487(5) 


0.4021(9) 




-0.310(2) 


1.340(2) 


1.7 


5 


1.224(1) 


0.8169(9) 




-0.177(9) 


1.473(9) 


1.9 


5 


1.375(2) 


0.727(1) 




-0.20(1) 


1.45(1) 


2.2 


5 


1.761(2) 


0.5680(5) 




-0.294(4) 


1.356(4) 


2.4 


5 


2.180(6) 


0.459(1) 




-0.308(4) 


1.342(4) 


TABLE IV. Results from fits to the p = 


1.7, ^ = 3 g 


iueball correlators. N^xp is 


the number of 


exponentials used in each fit, and the fit ranj 


je refers to the temporal separations tmin — tmax used 


in each fit. Q is the standard quality of fit. 


Energies are given in 


^. Final mass estimates are 


highlighted in boldface. 












Channel 


^exp 


fit range 




Q 


energies 


overlaps 


At-" 


1 


0-5 


0.71 


0.58 


1.061(7) 


1.000(4) 




1 


1-5 


0.83 


0.48 


1.05(2) 


0.99(2) 




1 


2-5 


1.20 


0.30 


1.04(5) 


0.96(9) 


E++ 


1 


0-3 


0.02 


0.98 


1.653(8) 


1.000(3) 




1 


1-4 


0.12 


0.88 


1.65(5) 


1.00(5) 




2 


1-3 


1.37 


0.25 


1.65(5) 


0.98(5) 












2.5(2) 


1.1(2) 




1 


0-3 


0.60 


0.55 


1.99(1) 


1.000(2) 




1 


1-3 


0.06 


0.81 


1.92(6) 


0.94(6) 




2 


1-3 


0.65 


0.58 


1.90(6) 


0.88(8) 












2.5(3) 


0.8(3) 


Tr 


1 


1 - 3 


0.62 


0.43 


2.09(9) 


0.95(8) 
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TABLE V. Results from fits to the /? = 1.9, ^ = 3 glueball correlators (see Table 



Channel 




fit range 


xVdof 


Q 


energies 


overlaps 




1 
1 

2 


0- 4 

1- 4 
1 - 3 


1.05 
1.46 
0.74 


0.37 
0.23 
0.53 


0.878(6) 
0.87(1) 

0.88(1) 
1.8(1) 


1.000(5) 
0.99(1) 
1.00(1) 
1.0(1) 




1 
1 


0- 4 

1- 4 

1 Q 

1 — 3 


0.11 
0.04 
O.Ub 


0.95 
0.96 

n no 


1.493(9) 
1.47(4) 

1.9(1) 


1.000(3) 
0.98(3) 
0.9o(oj 
0.91(9) 




1 
1 

2 


0-3 
1 - 3 
1 - 3 


1.74 
0.22 
0.23 


0.18 

0.64 
0.88 


1.681(8) 
1.60(4) 

1.58(4) 
2.6(2) 


1.000(2) 
0.92(4) 
0.88(4) 
1 3f3l 


r+- 


1 


1 - 3 


0.24 


0.63 


1.80(6) 


0.90(5) 


LAULtj Vi. 


Results from fits to the (3 


= 2.0, ? = 3 


glueball correlators (see Table [ 


LYI). 


Channel 


^exp 


fit range 


xVdof 


Q 


energies 


overlaps 




1 
1 

2 


0-6 
1 - 6 
1 - 3 


0.94 
1.15 
2.33 


0.45 
0.33 
0.07 


0.794(4) 
0.797(9) 

0.794(8) 
1.61(5) 


1.000(4) 
1.002(9) 
1.001(8) 
0.95(5) 




1 
1 
2 


0- 4 

1- 4 
1 - 3 


0.67 
0.42 

1 .ou 


0.57 
0.65 
0.27 


1.423(6) 
1.40(2) 

1.39(2) 
1.81(6) 


1.000(2) 
0.97(2) 

U.c/Ul ^ I 

0.94(5) 


rn++ 


1 
1 
2 


0- 4 

1- 4 
1 - 3 


2.54 
0.67 
0.53 


0.05 
0.51 
0.66 


1.559(5) 
1.49(2) 

1.49(2) 
2.11(7) 


1.000(2) 
0.94(2) 
0.94(2) 
1.01(7) 




1 


1-4 


0.02 


0.98 


1.68(3) 


0.92(3) 
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TABLE VII. Results from fits to the P = 2.2, ^ = 3 glueball correlators (see Table 



Channel 


Nexp 


fit range 


xVdof 


Q 


energies 


overlaps 




1 


0-7 


1.12 


0.35 


0.659(4) 


0.998(6) 




1 


1-7 


0.76 


0.58 


0.649(8) 


0.988(8) 




2 


1-4 


1.02 


0.41 


0.647(8) 
1.25(3) 


0.984(8) 
0.93(2) 




1 


1-4 


0.10 


0.90 


1.19(2) 


0.95(2) 




1 


2-4 


0.10 


0.75 


1.17(6) 


0.9(1) 




2 


1 - 3 


0.32 


0.81 


1.19(2) 
1.62(4) 


0.95(2) 
0.99(4) 




1 


0-4 


2.32 


0.07 


1.280(6) 


1.001(2) 




1 


1-4 


0.92 


0.40 


1.24(2) 


0.96(2) 




1 


2-4 


0.02 


0.88 


1.16(6) 


0.8(1) 




2 


1 - 3 


1.30 


0.27 


1.24(2) 
1.65(5) 


0.96(2) 
0.87(4) 




1 


1 - 4 


1.13 


0.32 


1.48(3) 


0.98(3) 



TABLE VIII. Results from fits to the /3 = 2.4, ^ = 3 gluebah correlators (see Table 



Channel 




fit range 


xVdof 


Q 


energies 


overlaps 




1 


1 - 8 


0.19 


0.98 


0.548(6) 


0.988(6) 




1 


2-8 


0.22 


0.96 


0.550(9) 


0.99(2) 




2 


1-4 


1.12 


0.35 


0.550(6) 


0.991(6) 












1.03(2) 


0.96(1) 


E++ 


1 


0-5 


1.30 


0.27 


1.012(4) 


1.000(3) 




1 


1 - 5 


0.31 


0.82 


0.995(9) 


0.982(9) 




2 


1 - 3 


2.14 


0.09 


0.993(8) 


0.982(8) 












1.37(2) 


1.00(2) 


-'2 


1 


0-5 


3.98 


0.00 


1.035(4) 


1.001(2) 




1 


1 - 5 


0.28 


0.84 


1.006(8) 


0.969(9) 




2 


1 - 3 


0.53 


0.66 


1.006(8) 


0.966(8) 












1.42(2) 


0.99(2) 




1 


1-4 


0.75 


0.47 


1.24(1) 


0.98(1) 
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TABLE IX. Results from fits to the /3 = 2.6, ^ = 3 glueball correlators (see Table |^). 



Channel 




fit range 




Q 


energies 


overlaps 




1 
1 

2 


1 - 10 
2-10 
1-4 


0.88 
1.01 
0.89 


0.53 
0.42 
0.50 


0.464(7) 

0.46(1) 

0.464(6) 

0.84(1) 


0.986(8) 
0.98(2) 
0.988(8) 
0.96(1) 




1 
1 


1 - 6 
2-6 

i — 4 


0.60 
0.36 

n on 
0.89 


0.67 
0.78 
O.oU 


0.781(8) 

0.76(2) 
0. Iii2[\)) 
1.09(2) 


0.984(8) 
0.94(4) 

0.9oz(9j 

0.95(2) 




1 
1 

2 


1 - 6 
2-6 
1-4 


0.31 

0.12 
0.53 


0.87 

0.95 
0.79 


0.777(8) 

0.76(2) 

0.777(8) 

1.12(1) 


0.977(8) 
0.95(3) 
0.976(8) 
0.99(1) 


r+- 


1 


1 - 5 


0.27 


0.85 


0.97(1) 


0.97(1) 


iAJDLrv A. 


Results from fits to the [3 


= 1-7, ? = 5 


glueball correlators (see Table [ 




Channel 


^exp 


fit range 


xVdof 


Q 


energies 


overlaps 




1 
1 

2 


0- 5 
1 - 5 

1- 4 


1.36 
0.26 
0.38 


0.24 
0.86 
0.89 


0.585(3) 
0.578(5) 

0.578(5) 
1.19(2) 


0.999(4) 
0.992(5) 
0.992(5) 
0.97(2) 




1 
1 
2 


0- 5 
1 - 5 

1- 4 


2.09 
0.12 
0.29 


0.08 
0.95 
0.94 


0.943(3) 
0.924(8) 

0.924(8) 
1.29(2) 


1.000(2) 
0.981(7) 
0.979(7) 
0.98(1) 




1 
1 
2 


0- 5 
1 - 5 

1- 4 


1.25 

i.Oo 

0.87 


0.29 
0.19 
0.52 


1.107(3) 
1.103(8) 

1.104(9) 
1.41(2) 


1.001(2) 
u.yy / yij ) 
0.997(9) 
0.94(2) 




1 
1 
2 


0-3 
1 - 3 
1 - 3 


2.65 
0.21 
0.39 


0.07 
0.65 
0.76 


1.214(4) 
1.19(1) 

1.18(1) 
1.55(3) 


1.000(2) 
0.97(1) 
0.97(1) 
0.92(2) 
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TABLE XI. Results from fits to the /3 = 1.9, ^ = 5 glueball correlators (see Table |^). 



Channel 


Nexp 


fit range 




Q 


energies 


overlaps 




1 


1 - 9 


1.26 


0.26 


0.475(4) 


0.992(5) 




1 


2-9 


1.11 


0.35 


0.468(6) 


0.98(1) 




2 


2-4 


1.03 


0.38 


0.468(6) 
0.92(3) 


0.98(1) 
0.85(5) 




1 


1 - 6 


0.42 


0.80 


0.844(6) 


0.992(6) 




1 


2-6 


0.20 


0.90 


0.83(1) 


0.97(2) 




2 


1-4 


0.69 


0.66 


0.844(6) 
1.09(1) 


0.992(6) 
0.95(1) 




1 


1 - 5 


0.91 


0.43 


0.918(7) 


0.982(6) 




1 


2-5 


1.11 


0.33 


0.91(2) 


0.96(3) 




2 


1-4 


1.44 


0.19 


0.918(6) 
1.23(1) 


0.981(6) 
0.96(1) 




1 


1 - 5 


0.30 


0.83 


1.053(8) 


0.979(7) 




1 


2-5 


0.18 


0.84 


1.04(2) 


0.95(4) 




2 


1-4 


0.48 


0.82 


1.052(9) 
1.30(1) 


0.97(1) 
0.92(1) 



TABLE XII. Results from fits to the /3 = 2.2, ^ = 5 glueball correlators (see Table IV). 



Channel 




fit range 


xVdof 


Q 


energies 


overlaps 




1 


1 - 14 


0.60 


0.84 


0.362(3) 


0.998(5) 




1 


2-14 


0.50 


0.90 


0.366(4) 


1.004(7) 




2 


1-4 


0.86 


0.52 


0.362(3) 
0.697(6) 


0.998(5) 
0.970(7) 


E++ 


1 


1-7 


1.33 


0.25 


0.667(4) 


0.982(4) 




1 


2-7 


1.65 


0.16 


0.666(7) 


0.98(1) 




2 


1-4 


0.57 


0.75 


0.667(4) 
0.878(7) 


0.980(5) 
0.968(7) 


r++ 


1 


1 - 8 


1.05 


0.39 


0.686(4) 


0.983(3) 




1 


2-8 


1.20 


0.31 


0.683(6) 


0.98(1) 




2 


1-4 


0.77 


0.59 


0.686(3) 
0.938(5) 


0.982(3) 
0.970(5) 




1 


1 - 6 


0.48 


0.75 


0.819(4) 


0.974(4) 




1 


2-6 


0.57 


0.63 


0.82(1) 


0.98(2) 




2 


1-4 


0.60 


0.73 


0.820(4) 
1.025(8) 


0.974(4) 
0.956(7) 
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TABLE XIII. Results from fits to the (3 = 2.4, = 5 glueball correlators (see Table 



Channel 


TVexp fit range 


xVdof 


Q 


energies 


overlaps 




1 


1 - 13 


1.35 


0.19 


0.303(3) 


0.995(7) 




1 


2-13 


1.15 


0.32 


0.307(4) 


1.000(8) 




2 


1 - 5 


0.86 


0.56 


0.304(3) 


0.994(7) 












0.569(4) 


0.972(5) 




1 


1 - 9 


1.33 


0.23 


0.538(3) 


0.992(3) 




1 


2-9 


1.46 


0.19 


0.536(5) 


0.986(9) 




2 


2-5 


2.71 


0.01 


0.538(5) 


0.99(1) 












0.723(9) 


0.94(1) 


rr+ + 


1 


1 - 9 


1.08 


0.38 


0.542(2) 


0.988(3) 




1 


2-7 


0.34 


0.85 


0.540(4) 


0.982(7) 




2 


2-5 


1.28 


0.26 


0.538(4) 


0.978(8) 












0.730(8) 


0.94(1) 




1 


2-6 


1.03 


0.38 


0.652(5) 


0.95(1) 




2 


z — D 


1.77 


0.07 


0.648(6) 


0.95(1) 












0.794(9) 


0.88(2) 


TABLE XIV. 


Summary of final 


mass estimates from all - 


= 3 and ^ = 5 simulations. 


P 




atm{At+) 


atm{E++) 




atm(r2++) 


atm{T^^) 


1.7 


3 


1.05(2) 


1.65(5) 




1.92(6) 


2.09(9) 


1.9 


3 


0.87(1) 


1.47(4) 




1.60(4) 


1.80(6) 


2.0 


3 


0.797(9) 


1.40(2) 




1.49(2) 


1.68(3) 


2.2 


3 


0.649(8) 


1.19(2) 




1.24(2) 


1.48(3) 


2.4 


3 


0.548(6) 


0.995(9) 




1.006(8) 


1.24(1) 


2.6 


3 


0.464(7) 


0.781(8) 




0.777(8) 


0.97(1) 


1.7 


5 


0.578(5) 


0.924(8) 




1.103(8) 


1.19(1) 


1.9 


5 


0.475(4) 


0.844(6) 




0.918(7) 


1.053(8) 


2.2 


5 


0.362(3) 


0.667(4) 




0.686(4) 


0.819(4) 


2.4 


5 


0.303(3) 


0.538(3) 




0.542(2) 


0.652(5) 
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TABLE XV. 
volumes. 


Glueball mass estimates 


in terms of a^ ^ 


for (3 = 2.4, i = 


3 and various lattice 


Channel 


Ls/as = 4 


Ls/as = 5 


Ls/as = 6 


Ls/as = 8 


^*++ 

171+ + 

T2++ 

2^2*++ 

ri+- 


0.483(6) 

0.863(9) 

0.848(9) 

1.11(2) 

1.032(9) 

1.37(2) 

1.245(13) 


0.551(6) 

0.97(1) 

0.980(10) 

1.27(2) 

1.005(9) 

1.42(2) 

1.235(11) 


0.545(8) 

1.00(3) 

1.004(19) 

1.41(4) 

1.002(14) 

1.40(3) 

1.209(25) 


0.548(6) 

1.027(17) 

0.995(9) 

1.37(2) 

1.006(8) 

1.42(2) 

1.237(14) 


TABLE XVI. The effects of simulating in a finite box: results from fits of Eq. 16 to the energy 
estimates given in Table XV from lattices of spatial extent Lg/as = 4,5,6,8. The fit also 
includes an energy estimate for Ls/as = 3. The final column estimates the expected finite-volume 
errors in glueball masses from the Ls/as = 8 simulation at /? = 2.4 and ^ = 3. These errors are 
estimated by mG{8uj^) / mcioo) — 1 using Eq. [l^. 


Channel 


atmcioo) 


Ag 


xVdof 


% correction 


E++ 

rp+ + 

n- 


0.554(4) 
1.002(7) 
1.003(6) 
1.223(7) 


260(37) 
319(23) 
-59(24) 
-66(4) 


2.3 
0.62 
0.15 
1.05 


-0.020 
-0.024 
+0.004 
+0.005 
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TABLE XVII. Extrapolations of the glueball mass estimates to the continuum hmit for the 
^ = 3 runs. The three scahng forms fo, ip2, and 994 which are fit to the data are given in Eqs. 17, 
18 and 19. The values indicated in bold are taken as our final continuum mass estimates. 



( , n 51 n Ti p>l 


Trl" fllTlpl'IOTl 


' 0'"'G 






/dof 












n 71 






Ql (7) 


— n 96i'9n'i 




n 49 




( n A 






—0 S3l'2.^'l 


46 




^0 














5.66(7) 


1 n2f2i'i 




0.90 




( A 






1 29r27'l 


0.87 


r+- 




7.22(5) 







1.81 




f2 


7.44(11) 


-0.71(31) 




0.96 




ip4 


7.32(7) 




-0.87(40) 


1.07 


E*++ 


^0 


7.99(8) 






2.44 




<P2 


8.52(20) 


-1.8(6) 




0.31 




994 


8.25(12) 




-2.6(9) 


0.12 


rn*++ 


^0 


8.37(6) 






1.33 




(P2 


8.32(16) 


0.2(7) 




1.94 




(P4 


8.34(8) 




0.5(11) 


1.89 
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TABLE XVIII. Extrapolations of the glueball mass estimates to the continuum hmit for the 
^ = 5 runs. The three scahng forms (po, ip2, and 934 which are fit to the data are given in Eqs. 17, 
18 and 19. The values indicated in bold are taken as our final continuum mass estimates. 



Channel 


fit function 




rouiG 


C2 




C4 


xVdof 




VO 
¥'2 

</34 




5.98(5) 
5.91(3) 


-0.41(12) 




-0.50(14) 


^ 03 
1.55 
0.81 




<4>2 

(/74 




fi 1 ^C?) 
5.53(4) 
5.82(3) 


1.69(11) 




2.02(13) 


89 
5.16 
1.32 


r+- 

1 






7.21(2) 

7.09(7) 
7.15(4) 


0.31(16) 




0.32(18) 


1.61 
0.57 
0.80 


E*++ 


^0 
V2 

(P4 




7.67(4) 

8.12(14) 

7.91(8) 


-1.2(4) 




-1.5(4) 


5.88 
0.02 
0.27 


rn*++ 


^0 
^2 
(P4 




8.29(3) 
7 Rnn 

( .oU\^iU ) 

8.06(6) 


1.3(3) 




1.4(3) 


9.17 
1.37 
2.48 


TABLE XIX. Extrapolations of 
fit functions used are given in Eqs. 
mass estimate for the scalar glueba 


the scalar 
20 and 21 
1. 


glueball mass estimates to the continuum limit. The 
. The value indicated in bold is our final continuum 


State 


fit function 


roniG 




C2 


C4 


CL 


xVdof 




<flL 


3.98(15) 

3.86(8) 


-18(4) 18(5) 
-3.5(4) 4.5(5) 


0.96(13) 


0.25 
0.55 




9^2,4 


6.93(19) 


-5.3(1.1) 8.^ 


^(1.4) 




0.41 
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TABLE XX. Estimates of using results from various quenched lattice simulations with the 
Wilson gluonic action. The simple average r^^ = 410(20) MeV of the last column is taken as our 
estimate. 



Source /? Quark Action Scale Setting Quantity a ^ (GeV) Tq ^ (MeV) 



NRQCD [H 


6.0 


NRQCD 


T{2S-1S, lP-15) 


2.4(1) 


434(23) 


NRQCD [H 


5.7 


NRQCD 


J/^P{1P-1S) 


1.23(4) 


430(16) 


LANL |26| 


6.0 


Wilson 


Mp 


2.330(41) 


422(16) 


GFU ^ 


6.17 


Wilson 




2.93(11) 


419(17) 


JLQCD p8| 


6.1 


Wilson 


J/i){lP-lS) 


2.54(7) 


394(13) 


JLQCD p8| 


6.3 


Wilson 


J/il){lP-lS) 


3.36(11) 


401(14) 


JLQCD p9\ 


6.3 


Wilson 


Mp 


3.41(20) 


406(24) 


FNAL |3C1 


6.1 


FNAL 


J/V'(1P-15) 


2.55(8) 


395(14) 


BLS |31j 


6.3 


heavy-light 




3.21(9) 


383(11) 


ELS 131] 


6.3 


heavy-light 


Mp 


3.44(9) 


410(11) 


UKQCD II 


6.2 


Wilson 


Mp 


2.77(16) 


376(22) 
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